Absence of a Phase Transition in a Three— Dimensional Vortex 

Glass Model with Screening 
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Abstract 

We study the gauge glass model for the vortex glass transition in type- 
II superconductors, including screening of the interaction between vortices. 
From the size dependence of the domain wall energy we find that, in two- 
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, dimensions, the transition is at T = both with and without screening but 
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Cp I the exponents are different in the two cases. In three-dimensions, we find 
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I that screening destroys the finite temperature transition found earlier when 

^ I screening was neglected. 
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The question of whether the hnear resistance, i.e. the resistance in the hmit of a van- 
ishingly small current is really zero in a type-II superconductor in fields greater than 
ifci has received considerable attention recently, both theoretical and experimental, because 
fluctuation effects are so much bigger for high Tc compounds than for low Tc materials. 
In fact much of the H — T phase diagram of high materials is occupied by a "vortex 
liquid" regime in which the resistance has dropped because superconducting short range 
order has formed, so flux lines exist locally, but the resistance is not yet zero because the 
flux lines move under the action of a Lorentz force due to the current, and hence give rise 
to a voltage 0. An important question, then, is whether, at lower-T, the flux lines can be 
collectively pinned by defects sufficiently strongly that they have no linear response to the 
Lorentz force, in which case the linear resistance is zero. Such as state was first proposed for 
bulk superconductors by M. P. A. Fisher ^ and called the vortex glass. Note that disorder, 
which is essential to get a vanishing resistance, destroys the flux lattice ^ predicted by 
mean field theory. 

It is important to distinguish whether the dominant type of disorder is point-like, which 
could be caused by oxygen vacancies in high-Tc materials, or whether it is due to extended 
defects such as twin boundaries or artificially created columnar defects, The latter 
case is analogous to bosons moving in a random potential in two dimensions 0. For this 
problem the vortex glass state is called a Bose glass [§], which undoubtedly occurs in 
three dimensions, though some controversies still remain about the details. However, the 
situation is less clear for point disorder, the case studied here. Much of the theoretical 
work has studied a simplified model, called the gauge glass. In two dimensions, this has 
a transition at T = [|10|,|ll|, whereas three dimensions seems close to the lower critical 
dimension ||12|,|lT|, though there does appears to be a finite |T^ in this case. Several 



experiments |T3| on samples of YBCO have found reasonably convincing evidence for a 
phase transition in non-zero field, presumably to some sort of vortex glass state. 

Studies of the vortex glass transition with point disorder have so far neglected screening 
of the (otherwise logarithmic) interaction between vortices. As we shall see below, screening 



cannot be neglected very close to Tc- The purpose of this paper is to see if the finite-Tc 
vortex glass transition survives the introduction of screening. We find, at least for the gauge 
glass model, that it does not. Hence, assuming that the gauge glass model is in the same 
universality class as experimental systems with point disorder, the observed "transition" in 
these systems would actually be rounded out sufficiently close to T^.. 

Screening will alter the critical behavior when the penetration depth. A, becomes equal 
to the correlation length, ^, because the important fluctuations on scale ^ will no longer be 
screened. The behavior of A and ^ are sketched in Figure 1. In the mean field regime, both A 
and ^ diverge as where t = T — so the ratio A/^ = k (> 1) is roughly constant and 

is a property of the material. As the temperature is reduced there will be a crossover from 
mean field behavior to the critical behavior of an uncharged superfiuid, in which screening 
can be neglected, at the Ginzburg temperature, 0. In this critical regime, ^ increases 
faster than A 0, so the two lengths eventually become equal at T = Tgc, say. Screening 
effects will therefore be important for < T < Tg^- Since C,{Tsc) ~ A(Tsc) ~ A^-, we have 
^{Tsc)/^{Tg) ~ 1^"^, see Figure 1. Now k, ~ 100 for high T^. compounds so the region where 
screening is important is small and has probably not been accessed in experiments so far. 
Nonetheless we feel that it is useful to study the behavior in this region, because (i) this 
regime may be accessible in experiments on certain materials and ranges of field, and (ii) 
theoretically, it is interesting to know whether or not the linear resistance is really zero. 

Next let us discuss the gauge glass which we model used to study effects of screening. 
Including a fluctuating vector potential, the Hamiltonian is 



The phase of the condensate, (f)i, is deflned on each site, /, of a regular lattice, square for 
two dimensions, and simple cubic for d = 3, with N = L'^ sites. For now we assume periodic 
boundary conditions. The sum is over all nearest neighbor pairs on the lattice. The effects of 
the external magnetic fleld and disorder are represented by the quenched vector potentials, 
Ajj, which we take to be independent random variables with a uniform distribution between 




(1) 
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and 27r. The fluctuating vector potential on each link, Aij, is integrated over from — oo to 



oo . In the last term, which is the usual magnetic energy, the sum is over all elementary 
squares on the lattice, and the curl, V x ^, is the directed sum of the vector potentials 
round the square, taking into account that Aij = —Aji. The "charge", e, is a measure of 
the coupling between the condensate and the fluctuating vector potential. The previously 
studied model |jl0|,|l2|, which neglects screening, corresponds to e = 0. 

It is technically easier to study this model in the vortex representation rather than the 
above phase representation. To do this we replace the cosine by the periodic Gaussian 



(Villain) function ||T5l, and perform standard manipulations [|15],|16|. For our purposes it is 
also essential to incorporate the periodic boundary conditions, which leads to an additional 



term |T0| , p!7| . For two dimensions, we find that the vortex Hamiltonian is 

nv = -^E(^^ - h)G{^-J){n, - b,) + E \/($" - 6"^P^) , (2) 

i,j a=x,y 

where the {n,} run over all integer values, subject to the "charge neutrality" constraint 
J2i^i = 0, and are interpreted as the vortices, G{i — j) is the screened vortex interaction. 



G{i - j) ^ ]_Y- l-exp[zk- (r; - rj)] 
(27r)2 4 - 2 cos A;^. - 2 cos fcy + e2 



(3) 



and V{u) is the Villain potential given by 

'-V{u " 



exp 
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J2 exp 



m=—oo 



(4) 



V{u) is periodic with period 27r, and, at T = 0, is equal to v? /2 for — tt < m < tt. The 
vortices sit on the sites, z, of the dual lattice which are in the centers of the squares of the 
original lattice. The hi are given by (l/27r) times the directed sum of the quenched vector 
potentials on the links of the original lattice which surround the site i of the dual lattice. 



As usual e^^ 



0. The function $^ is defined by 



(5) 



and similarly for Here we indicate explicitly the x and y coordinates of the lattice 
points: we use integers for the dual lattice, x = 1, 2, . . . , L, so the coordinates of the original 



lattice are half-integers. The link variable Arc{x + 1/2, 1/2) connects sites (x + 1/2, 1/2) and 
(x + 3/2, 1/2). The function P is 2tt/L times the "dipole moment" of the vortex and flux 
distribution, i.e. 

P = f Er Kr)-Kr)] ■ (6) 

The Hamiltonian is independent of the choice of origin for the coordinate system: choosing a 
new origin alters $ and P but it is straightforward to show that the arguments of the Villain 
functions only change by a multiple of 27r and so the functions themselves are unchanged. 
Notice that for e = the interaction between vortices is logarithmic but that for e 7^ this 
logarithmic variation is screened beyond a distance of order 1/e. 



In d = 3 the vortices are integer variables which lie on the links of the dual lattice [jT6 
We find that the vortex Hamiltonian, including periodic boundary conditions, is 

Hy = -VK-b,)-(n,-b,)G(z-j)+ Y: V - le^p.Cf'A , (7) 

where the rij are subject both to the global constraint J2i = 0, analogous to d = 2, and 
to the local constraint, [V • n]j = for each site |TB[. The vortex interaction, G{i — j), is 



obtained by replacing the RHS of Eq. by its three-dimensional analogue, C is given by 
= § E ^ N^r) - yq^ir)] + f E ^ [^l'ir)Sy,L - yq^{r)6.,L] , (8) 

where q(r) = n(r) — b(r), with similar expressions for the other components, and # is 
defined as in Eq. (|^). The physical interpretation of C is quite simple. Consider the change 
in C when 1 is added to all the rij on a closed loop in the xy plane. Eq. (||) tells us that C^^ 
is then increased by 271/1? times the area of the loop, while C^^ and C^^ are unchanged. 
It is easy to see that this is true if the loop does not go along links [|T^ at which there is 
a discontinuity in our labeling of the coordinates of the sites, because only the first term 
in Eq. then contributes, but it is not complicated to show that the result is true for 
all loops. A change in the boundary conditions, from, say, periodic to anti-periodic in in 
the z-direction, can be incorporated by adding vr/L to all the Az, and so has the effect of 



changing $^ by vr. This can be compensated for, in the arguments of the Villain functions, 
by forming a vortex loop in the xy plane whose area is or a combination of loops with 

the same total area. A derivation of Eqs. and (|^) will be given elsewhere . 



The gauge glass model has random fluxes with no preferred direction, as distinct from 
the net applied field which occurs in experimental systems. This difference is irrelevant in 
an expansion in 6 — e dimensions, and so may also be irrelevant in = 3, though a firm 
demonstration of this is lacking. 



To investigate whether a transition occurs we calculate the defect wall energy pO|. For 
a given realization of the disorder, we compute the difference in ground state energy, AE, 
between the system with periodic boundary conditions, and with boundary conditions which 
are anti-periodic in one direction and periodic in the others. To find the ground state, we 
start from different random initial configurations and quench down to the nearest local 
energy minimum. We repeat this many times, keeping track of the lowest energy found so 
far. We monitor the r.m.s. energy difference, AErms, discussed below, after a logarithmically 
increasing number of quenches, 1, 3, 10, . . . and stop when the last two estimates are virtually 
identical. For a number of cases we ran for much longer than this and verified that the results 
were unchanged. We also also performed some calculations in the phase representation, using 
the Villain form for the interaction, and checked that identical results were obtained. 

For frustrated systems, such as the model we study here, the average of AE over all 
samples will be zero, so we look at the root mean square average, AErms = [^-^^]av^ > where 
[. . .]av denotes an average over samples. We investigate the size dependence of AErms and 
define an exponent 6 by 

AErms -^L' . (9) 

If < then large domains cost little energy so, at any finite T, the system will break up 
into domains and long range order will be destroyed. By equating the domain wall energy 
on the scale of the correlation length with ksT one finds that the correlation length varies 
with T as with v = l/\9\. If, on the other hand, 6' > then then the rigidity of the 

6 



ordered state should persist to finite temperature. 

Our results for d = 2 are presented in Figure 2 for 2 < L < 8. It was not possible to get 
convergent results for AErms for larger sizes. For e = the slope is ^ ~ —0.5, which is in 



good agreement with earlier work |T^, and gives a correlation length diverging as roughly 
Experimental evidence for this has recently been obtained on very thin films of 
YBCO. For e 7^ the slope is clearly more negative, indicating that screening is a relevant 
perturbation. 

The curvature of the data indicates that the intermediate values of e are in a crossover 
regime. It is therefore be useful to study directly the large e limit. To do this, notice 
from Eq. that the vortex energy becomes very small in this limit, so the Villain term 
in Eq. (Q) dominates and acts as a constraint fixing the dipole moments of the vortex 
distribution to the values which minimize the Villain term. Also, in this limit, G{0) = 0, 
and G{i — j) = (27r/e)^ plus exponentially small corrections which we neglect. Because 
of charge neutrality, one can subtract the constant (27r/e)^ from all the interactions, with 
the result that only G{0) (= — (27r/e)^) is non-zero. Hence, for e — 00, AErms^^ tends to 
a finite value which can be found by by neglecting the Villain terms in the energy, setting 
G{0) = — (27r)^ and G{i—j) = for i 7^ j, and performing a constrained minimization fixing 
not only the total "charge" but also the total dipole moment of the vortex distribution. On 
changing the boundary conditions, the x-component of the dipole moment is changed by 
L/2 [^]. Results for e'^AE^ms obtained in this way are shown in the inset to Figure 2. The 
data follows a straight line on a log-log plot with a slope of ^ —1.36, which is also roughly 
the slope of the data for the larger values of e and L in the main part of the Figure. 

This result is surprising because there are simple arguments which predict that 



6 = —2 in the large e limit. One way of seeing this is that, since the vortices don't interact, 
the probability that a vortex is excited out of its ground state is obtained just by comparing 
its on-site excitation energy with ksT. Since there is a finite density of states, the density of 
excited vortices is of order /c^T, and hence a typical separation (which we interpret to be the 
correlation length) of order T~^/^, and so 9 = = 2. Assuming that this is correct, the 



different result found here must reflect the importance that the constraints (total charge = 0, 
fixed dipole moment) have for the rather small system sizes studied here. Note, though, that 
the presumably correct result, 9 = —2, predicts an even /aster drop of AErms than found in 
the numerics, and that the numerics is correct in predicting that the screened system is in 
a different universality class with a weaker divergence at T = 0. 

Our results for AErms in d = 3 are presented in Figure 3. The range of sizes, L = 2,3 
and 4, is very small, but analogous calculations on similarly small sizes provided one of 
the first clear indications for the, apparently correct, result that Tc is finite in 3-d Ising 



spin glasses ||20|. For e = 0, AErms is roughly independent of size, in agreement with earlier 



work |T2|, which indicates that the system is close to the lower critical dimension. For e 7^ 0, 
however, AErms, decreases with increasing size, this effect being more pronounced for larger 
values of e. One can also study the large e limit, as for d = 2, but we were only able to 
get convergence for L = 2 and 3, and, with only two points, one cannot check for straight 
line behavior of the log-log plot. Although the sizes are very small, the results of Figure 3 
clearly imply that there is no finite temperature transition in d = 3. The finite temperature 
transition, for which some evidence was presented in Ref. []T2[, is therefore rounded out when 



one includes screening. 

To conclude, the presumed finite Tc vortex glass transition in the three-dimensional gauge 
glass model is rounded out by screening effects. The same will be true in real experimental 
systems with point disorder, assuming that the model is in the same universality class, which 
is plausible but not yet firmly established. Strictly speaking, the linear resistance would then 
not vanish, though it would become extremely small because the rounding only occurs very 
close to Tc. It would obviously be interesting to investigate this model by finite temperature 
Monte Carlo simulations, where larger sizes can be studied, and also to simulate a more 
realistic model with a net field in a particular direction. 

We have benefited enormously from stimulating discussions with M. P. A. Fisher. We 
would also like to thank D. A. Huse for valuable comments. This work is supported by NSF 
grants DMR 91-11576 and 94-11964. 
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FIGURES 

FIG. 1. A sketch of the log of the correlation length, ^, and penetration depth, A, against 
logi, where t = T — Tc, for a type-II superconductor. The slopes of the curves are indicated: is 
the correlation length exponent in the part of the critical region where screening can be neglected. 
Two crossovers occur. The first, at the reduced Ginzburg temperature, tc = Tq — Tc is when 
critical fluctuations first start to be important. The second, at tgc = Tg^ — ?c, is when A and ^ 
become comparable. Screening effects are important for t < tgc- For T > Tc A/^ = k (> 1), is 
constant, while for T^c < T < Tq, A^/^ is constant Q. Between Tq and T^c, increases by a factor 
of and A by a factor of n. 

FIG. 2. A log-log plot of the r.m.s. domain wall energy, against L for d = 2 for different 
values of the coupling e. The number of samples varied between 4000 for the smallest sizes to 200 
for the largest sizes and values of e. For e = the slope, 0, is approximately —0.5 in agreement 



with earlier work [10|. For e 7^ the curves bend down more steeply indicating that screening is 
a relevant perturbation. The insert shows a log-log plot of e^/S.E against L for d = 2 in the limit 
e — > 00. The number of samples varied between 2000 for L = 2 and 500 for L = 6. The line is a 
least squares fit and has a slope of —1.36. 

FIG. 3. A log-log plot of /S.E against L for d = 3 for different values of the coupling e. The 
number of samples was 4000 for L = 2 and at least 1000 for L = 4. For e = the slope is close to 



zero, in agreement with earlier work |12]. For e 7^ the curves bend down indicating that screening 



is a relevant perturbation which destroys the (possible) finite temperature transition at e = 0. 
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